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ABSTRACT 

We present a new, objectively defined catalog of candidate galaxy clusters based 
on the galaxy catalogs from the Digitized Second Palomar Observatory Sky Survey 
(DPOSS). This cluster catalog, derived from the best calibrated plates in the high lati- 
tude (|6| > 30°) Northern Galactic Cap region, covers 5,800 square degrees, and contains 
8,155 candidate clusters. A simple adaptive kernel density mapping technique, combined 
with the SExtractor object detection algorithm, is used to detect galaxy overdensities, 
which we identify as clusters. Simulations of the background galaxy distribution and 
clusters of varying richnesses and redshifts allow us to optimize detection parameters, 
and measure the completeness and contamination rates for our catalog. Cluster rich- 
nesses and photometric redshifts are measured, using integrated colors and magnitudes 
for each cluster. An extensive spectroscopic survey is used to confirm the photometric 
results. This catalog, with well-characterized sample properties, provides a sound basis 
for future studies of cluster physics and large scale structure. 

Subject headings: catalogues - surveys - galaxies: clusters: general - large-scale struc- 
ture of the Universe 

1. Introduction 

Over the past two decades, large numbers of astronomers have expended a great deal of effort on 
the study of galaxy clusters. Motivations for these works cover a wide range of astrophysical topics, 
from the study of galaxy evolution in dense environments, through the search for dark matter, the 
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measurement of the cluster mass function and its comparison to theoretical predictions, and on to 

the characterization of large scale structure in the universe. Despite this labor, understanding the 
underlying physics has often proved enigmatic, due to the complex physical nature of clusters, which 
is only compounded by the lack of an objectively selected, well-characterized statistical sample. 

To help remedy this situation, and provide a basis for future studies, we have undertaken the 
creation of a new cluster catalog which fulfills a number of fundamental criteria: 

1. Cluster detection should be performed by an objective, automated algorithm to minimize 
human biases and fatigue. 

2. The algorithm utilized should impose minimal constraints on the physical properties of the 
clusters, to avoid selection biases. 

3. The sample selection function must be well-understood, in terms of both completeness and 
contamination, as a function of both redshift and richness. The effects of varying the cluster 
model on the determination of these functions must also be known. 

4. The catalog should provide basic physical properties for all the detected clusters, such that 
specific subsamples can be selected for future study. 

In this paper, we describe how we generate a cluster sample that meets the above four criteria, 
and provide the final catalog for 5800n°. We discuss past cluster surveys and their limitations in 
§2. In §3 we briefly review the DPOSS survey (Djorgovski et al. 2003) which provides the basis for 
our catalog, while §4 describes the cluster detection technique (a modified version of that presented 
in Gal et al. (2000), hereafter Paper I). §5 describes simulations used both to optimize the detection 
algorithm, and define its statistical properties, including selection functions in richness and redshift 
for each plate. The results of these simulations are compared to extensive spectroscopic follow- 
up observations. Our photometric redshift estimator and richness measure are described in §6. 
An overview of the final catalog is provided in §7. We conclude with a brief discussion of future 
extensions of this survey (to cover the South Galactic Cap area, and the remaining areas of the 
NGP), and discuss some projects we have undertaken utilizing this sample. We use a cosmology 
with Hq = 67 km s^^ Mpc^^ and go = 0.5 throughout; other choices for these parameters have no 
significant effects on our results except for the poorest and most distant clusters. 



2. A Brief History of Cluster Surveys 

Surveys for galaxy clusters have only recently benefited from the automation afforded by 
computers. The earliest surveys relied on visual inspection of vast numbers of photographic plates, 
usually by a single astronomer. As early as 1938, Zwicky discussed such a survey based on plates 
from the 18" Schmidt telescope at Palomar, and results were presented in 1942 by both Zwicky 
and Katz &; Mulders. Even then, Zwicky, with his typical prescience, noted that elliptical galaxies 
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are much more strongly clustered than late-type galaxies. However, the true pioneering work in 

this field did not come until 1957, upon the publication of a catalog of galaxy clusters produced by 
George Abell as his Caltech Ph.D. thesis, which appeared in the literature the following year (Abell 
1958). Zwicky and collaborators followed suit a decade later, with their voluminous Catalogue of 
Galaxies and of Clusters of Galaxies (Zwicky, Herzog & Wild 1968). 

Abell used the red plates of the first National Geographic-Palomar Observatory Sky Survey. In 
selecting clusters, Abell applied a number of criteria in an attempt to produce a fairly homogeneous 
catalog. He required a minimum number of galaxies within two magnitudes of the third brightest 
galaxy in a cluster, a fixed physical area within which galaxies were to be counted, a maximum 
and minimum distance (redshift) to the clusters, and a minimum galactic latitude to avoid ob- 
scuration by interstellar dust. The resulting catalog, consisting of 1,682 clusters in the statistical 
sample, remained the predominant such resource until 1989. In that year, Abell, Corwin & Olowin 
(hereafter AGO) published an improved and expanded catalog, now including the Southern sky. 
These catalogs have been the foundation for many cosmological studies over the last four decades, 
despite serious questions about their reliability (which arc addressed later in this section). Some 
other catalogs based on plate material have also been produced, such as Shectman (1985), from 
the galaxy counts of Shane &: Wirtanen (1954), and a search for more distant clusters carried out 
on plates from the Palomar 200" by Gunn, Hoessel & Oke (1986). 

Only in the past ten years has it become possible to utilize the objectivity of computational 
algorithms in the search for galaxy clusters. These more modern studies required that plates be 
digitized, so that the data are in machine readable form. Alternatively, the data had to be digital 
in origin, coming from GGD cameras. Unfortunately, this latter option provided only small area 
coverage, so the hybrid technology of digitized plate surveys blossomed into a cottage industry, with 
numerous catalogs being produced in the past decade. All such catalogs relied on two fundamental 
data sets: the Southern Sky Survey plates, scanned with the Automatic Plate Measuring (APM) 
machine (Maddox, Efstathiou & Sutherland 1990) or GOSMOS scanner (to produce the Edin- 
burgh/Durham Southern Galaxy Catalog / EDSGC, Heydon-Dumbleton, Collins &; MacGillivray 
1989), and the POSS-I, scanned by the APS group (Pennington et al. 1993). The first objective 
catalog produced was the Edinburgh/Durham Cluster Catalog (EDGC, Lumsden, Nichol, Collins 
Sz Guzzo 1992), which covered 0.5 steradians ('^ 1,600 square degrees) around the South Galactic 
Pole (SGP). Later, the APM cluster catalog was created by applying Abell-like criteria to select 
overdensities from the galaxy catalogs, and is discussed in detail in Dalton, Maddox, Sutherland 
& Efstathiou (1997). The work by Odewahn & Aldering (1995), based on the POSS-I, provided 
a Northern sky example of such a catalog, while utilizing additional information (namely galaxy 
morphology). Some initial work on this problem, using higher quality POSS-II data, was performed 
by Picard (1991) in his thesis. 

In addition to these hybrid photo-digital surveys, smaller areas, to much higher redshift, have 
been covered by numerous deep CCD imaging surveys. Notable examples include the Palomar 
Distant Cluster Survey (PDCS, Postman et al. 1996), the ESO Imaging Survey (EIS, Olsen et 
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al. 1999, Lobo et al. 2000), Gonzalez, Zaritsky, Dalcanton, & Nelson (2001), KPNO/Deeprange 
(Postman, Lauer, Ocgerle, & Donahue 2002), and many others. None of these surveys provide 
the angular coverage necessary for large-scale structure and cosmology studies, and are typically 
designed to find rich clusters at high redshift. Only the Sloan Digital Sky Survey (SDSS, York et 
al. 2000) will provide large area, moderately deep CCD coverage. Cluster surveys from the SDSS, 
including those of Annis et al. (1999), Kim et al. (2002) and Goto et al. (2002), are only now 
appearing, and will cover a much smaller area until the photometric survey is completed. 

Despite these efforts, one thing is still missing: a catalog of galaxy clusters, produced by 
objective means, that is at least as deep as the Southern surveys, but which covers the Northern 
sky. This paper is intended to provide such a catalog. 

We note that there have also been many surveys for galaxy clusters at other wavelengths, 
most notably in the X-ray. All-sky surveys, such as RASS (Collins et al. 2000; Ebeling et al. 2000; 
Bohringcr et al. 2000) and EMSS (Gioia &: Luppino 1994), as well as pointed ROSAT observa- 
tions (Scharf et al. 1997; Romer et al. 2000), have been used to produce cluster catalogs at X-ray 
wavelengths. Future Sunyaev- Zel'dovich effect (SZ) surveys will also generate extremely important 
datasets. However, the relationship between optically and X-ray selected clusters is still not fully 
understood, so an optically selected catalog is essential. 

2.1. Limitations of Existing Catalogs 

Most of the optical studies to date have been limited by the statistical quality of the avail- 
able cluster samples. For instance, the Abell catalog suffers from a non-objective selection process, 
poorer plate material, a bias toward centrally-concentrated clusters (especially those with cD galax- 
ies), a relatively low redshift cutoff {z ~ 0.15; Bahcall & Soneira 1983), and strong plate-to-plate 
sensitivity variations. Still, many far-reaching cosmological conclusions have been drawn from it 
{i.e. Bahcall &; West 1992), although later studies have sometimes shown these to be flawed. 
Photometric errors and other inhomogeneities in the Abell catalog (Sutherland 1988; Efstathiou et 
al. 1992), as well as projection effects (Lucey 1983; Katgert et al. 1996) are serious and difficult- 
to-quantify issues. These effects have resulted in discrepant results on the correlation function 
(Bahcall & Soneira 1983; Dalton, Maddox, Sutherland & Efstathiou 1997; Miller et al. 1999), and 
later attempts to disentangle these issues relied on models to decontaminate the catalog (Suther- 
land 1988; Olivier et al. 1990). Also, Sutherland &; Efstathiou (1991) find a discrepancy between 
the angular and spatial correlation function of the Abell catalog, which is not found in the APM 
catalog (Dalton, Maddox, Sutherland & Efstathiou 1997). The extent of these effects is also sur- 
prisingly unknown; measures of completeness and contamination in the Abell catalog disagree by 
factors of a few. For instance. Miller et al. (1999) claim that under- or over-estimation of richness 
is not a significant problem, whereas van Haarlem, Prenk &; White (1997) suggest that one-third of 
Abell clusters have incorrect richnesses, and that one-third of rich {R > 1) clusters are missed. The 
largest study of i? > 1 clusters (Miller et al. 1999) suggests projection effects are not of great 
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concern for the Abell catalog; however, completeness cannot be gauged without deeper samples. 

Unfortunately, some of these problems will plague any optically selected cluster sample, including 
our own, but objective selection criteria and a strong statistical understanding of the catalog can 
mitigate their effects. 

Even some of the other objective catalogs preceding ours have their drawbacks. The APM 
group, for instance, used digitized J (blue) plates from the Southern Sky Survey; the use of a 
single blue band provides no color information to distinguish galaxy types, and is a poor choice 
for cluster detection because clusters are better delineated by redder, early-type galaxies in the 
redshift range probed {z < 0.3). Comparable surveys, such as the EDCC, already find a factor 
of two higher space density of clusters than Abell, and more sensitive CCD surveys find as many 
as five times more, although these results may be due to strong detection efficiency differences at 
lower richnesses. Other more recent surveys, such as the EDCCII (Bramel, Nichol & Pope 2000) 
have not yet achieved the area coverage of DPOSS. Additionally, the survey presented here utilizes 
at least one color (two filters) for photometric redshifts, and a significantly increased amount of 
CCD photometric calibration data. 



3. DPOSS: A Brief Overview 

The POSS-II (Reid et al. 1991) covers the entire northern sky {5 > —3°) with 897 overlapping 
fields (each 6.5° square, with 5° spacings), and, unlike the old POSS-I, has no gaps in the coverage. 
Approximately half of the survey area is covered at least twice in each band, due to plate overlaps. 
Plates are taken in three bands: blue-green. Ilia- J + GG395, Agg ~ 480nm; red, IIIa-F + RG610, 
AefT ~ 650nm; and very near-IR, W-N + RG9, XgS ~ 850nm. Typical limiting magnitudes reached 
are gj ~ 21.5, rp ~ 21.0, and in ~ 20.3, i.e., ~ 1"* - 1.5"* deeper than POSS-I. The image quality 
is improved relative to POSS-I, and is comparable to the Southern photographic sky surveys. 

The original survey plates are digitized at STScI, using modified PDS scanners (Lasker et 
al. 1996). The plates are scanned with 15/x (1.0") pixels, in rasters of 23,040 square, giving ~ 
1 GB/plate, or ~ 3 TB of pixel data total for the entire digital survey (DPOSS). Currently, 
astrometric solutions, provided by STScI, are good to rms ^ 0.5". 

An extensive effort to process, calibrate, and catalog the scans, with the detection of all objects 
down to the survey limit, and star/galaxy classifications accurate to 90% or better down to 1™ 
above the detection limit, has been undertaken over the past several years. Object detection and 
photometry is performed by SKICAT, a software system developed for this purpose (Weir et al. 
1995a,b,c), incorporating standard astronomical image processing packages, a commercial Sybase 
database, as well as a number of artificial intelligence and machine learning modules. Using this 
code, we measure ^ 60 attributes per object on each plate. Nearly all plates at |6| > 10° have 
been processed into catalogs. The catalogs are photometrically calibrated using extensive CCD 
sequences, with typical rms magnitude errors of 0.25™ at = 19.5; see Gal et al. (2003) for 
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details of the calibration procedure. A fraction of this error is due to systematic offsets in the 

photometric zero points between plates. Gal et al. (2003) show that the mean zero-point error is 
negligible, but has a la scatter of 0.07™ in the r-band, which can produce significant plate-to-plate 
depth variance. Our solution to this potential problem is discussed in §4. Star-galaxy separation 
is performed using a combination of FOCAS, neural networks, and decision trees, maintaining an 
accuracy of > 90% at < 19.5 (Odewahn et al. 2003). 

Each field in each band is processed individually. The three resulting catalogs are cross-matched 
to create a composite list of objects for the field. We require a detection in both the J and F bands 
so that we can measure the g — r colors of our galaxies; this also reduces the likelihood that any 
object is a spurious detection. Areas on the plate containing saturated objects are masked. These 
areas often contain large numbers of falsely identified galaxies, as the plate processing software, 
tuned to find faint objects, handles large, bright objects improperly. 

In this paper we utilize a total of 237 plates with good calibration (> 1, 000 calibrating galaxies 
per plate from the CCD fields), at 5 > 30°, in the North Galactic Cap region. The distribution of 
DPOSS fields (as well as the detected clusters) can be seen in Figure 14. The total area coverage 
is 5800 square degrees. In addition, we have run our procedures on one field (475) from the SGP 
region of DPOSS, which is used to obtain additional spectroscopic follow-up (see §5.3). 

4. The Detection Algorithm 

For this survey we use a modified version of the detection strategy first described in Paper I. 
We urge the reader to review Paper I for details of the adaptive kernel technique; we do not repeat 
the details here. In that earlier work, a color selection was first applied to the galaxy catalogs, after 
which we used an adaptive kernel technique (Silverman 1986) to produce galaxy density maps, and a 
bootstrap technique was used to create significance maps. The FOCAS object detection algorithm 
(Jarvis &; Tyson 1981) was then used to detect density enhancements in the significance maps, 
which we identified as candidate clusters. We generated cluster catalogs for only two sky survey 
fields, for which we performed extensive follow-up imaging and spectroscopy. Since the publication 
of that work, we have gained a greater understanding of the platc-to-plate variance in DPOSS, 
significantly modified our photometric calibration techniques, and learned a great deal about our 
detection algorithms as a result of the additional data obtained. This has led us to significantly 
modify our procedures, while maintaining the underlying principles and methodology. 

In Paper I, we utilized galaxies with magnitudes < 20.0. We have found that the ran- 
dom photometric errors, plate-to-plate zero-point offsets, and classification accuracy, can all reach 
unacceptable levels at that magnitude limit. Although some plates perform well to this limit, we 
found that imposing a slightly brighter magnitude limit, < 19.5, and requiring a g detection, 
produces a significantly more uniform galaxy catalog, without greatly sacrificing depth. Our intent 
with the current survey is to produce a catalog with good uniformity; therefore, we have elected to 
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be rather conservative in the data utilized. Future work will use fainter objects to create a higher 
redshift cluster catalog which may not be suitable for large scale structure work, but which will 
provide a useful sample for cluster studies over a larger distance/time baseline. Additionally, in 
Paper I we imposed a liberal selection in g — r color; we no longer apply this cut for the same 
reasons that we adopt the shallower magnitude limit. Because photometry in the g and r bands is 
completely independent, the use of a color cut effectively increases the pistoning due to zero-point 
offsets among plates by a factor -\/2. Finally, star/galaxy separation for DPOSS has been improved 
from the version used to produce the input data for Paper I, resulting in purer and more uniform 
galaxy samples being generated from the individual plate catalogs. The final catalog of objects 
now used in cluster detection therefore consists of all galaxies, with detections in both the g and 
r bands, having 15.0 < rur < 19.5. A typical galaxy catalog for a single plate contains ~ 27, 000 
galaxies over an area of ~ 34n°, resulting in a mean galaxy density of 0.79 x 10^ galaxies per square 
degree. 

The resulting galaxy catalog is used as input to the adaptive kernel (AK) density mapping 
algorithm. As described in Paper I, this technique uses a two-stage process to produce a density 
map. First, it produces an initial estimate of the galaxy density at each point in the map, which is 
then used to apply a smoothing kernel whose size changes as a function of the local density, with 
a smaller kernel at higher density. We refer the reader to Paper I for a more detailed description. 
As before, we generate our maps with 1' pixels. However, rather than the 900" kernel used earlier, 
we now use a significantly smaller kernel of 500" radius. Based on the simulations discussed in the 
next section, we found that this smaller kernel prevents over-smoothing the cores of higher-redshift 
{z ~ 0.3) clusters, while avoiding fragmentation of most low redshift {z ~ 0.08) clusters. The 
effect of varying the initial smoothing window is demonstrated in Figure 1. In this figure, we have 
placed four simulated clusters into a simulated background, representing the expected range of 
detectability in our survey. There are two clusters at low-z (0.08), and two at high-z (0.24), with 
one poor and one rich cluster at each redshift (100 and 333 total members, respectively). Of the 
100 (333) total cluster galaxies, only 54 (186) are brighter than our magnitude limit (m,. < 19.5) 
at z = 0.08, and 9 (17) at z = 0.24. The corresponding richnesses {Ngais, sec §6.2) are 80 and 25. 
The initial smoothing window is varied from 300" to 800" in 100" steps. This figure clearly shows 
the segmentation of rich, low- 2; clusters by small kernels, and the smoothing away of high-z clusters 
by large kernels. 

Unlike Paper I, we do not perform bootstrap tests to produce significance maps. Cluster 
candidates are detected directly from the AK map of the actual data. In the current work, we 
use SExtractor (Bertin & Arnouts 1996) to perform object detection (rather than FOCAS, as we 

did originally). We have found that SExtractor is both faster and more easily configured to meet 
our needs. For each plate, we find the set of SExtractor detection parameters (a pair of threshold 
level and minimum area) which produces an acceptable level of contamination by false clusters 
(10%). The simulations used to determine these parameters are discussed in §5.1. Once the set of 
optimal detection parameters for each plate is determined, further simulations are used to assess 
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the catalog completeness as a function of redshift and richness, as discussed in §5.2. We have opted 

to vary detection parameters between plates to maintain a constant level of contamination, even 
though the completeness may change. Due to the inhomogcncous nature of plate data, as well as 
large scale structure, setting a fixed threshold (in galaxy surface density) for all plates can lead 
to significant variance in both completeness and contamination between plates. By adjusting the 
parameters for each plate, we can minimize the variance of one of these quantities, and estimate 
the impact on the other. Bccrtisc contamination is more easily quantified (as an integrated rate, 
independent of richness and redshift) as the fraction of detections that are false, wc have chosen 
to adjust the parameters to maintain a fixed contamination rate. To allow large scale structure 
work, we provide a table (see Table 2) for each plate, including the RA/Dec boundaries used and 
the completeness function, as described in §5.2. 

Additionally, to assess the impact of photometric zero-point errors, and remove those cluster 
candidates whose detection is sensitive to these errors, wc generate a set of ten galaxy catalogs 
for each plate with random zero-point offsets added to the r-band magnitude, drawn from the 
known photometric error distribution for DPOSS given in Gal et al. (2003). Like the original data, 
these catalogs are clipped at 15.0 < rrir < 19.5, and an AK map generated for each one. Cluster 
detection is then performed using the same parameters as for the original map, and the resulting 
lists of cluster candidates are compared. Using the photometric redshifts (described in §6) for each 
candidate from the original map, we match clusters within a SOOkpc projected radius. In Table 
1 we present the fraction of total candidates that appear in a given number of maps. Only those 
candidates which appear in the original catalog, and more than 7 of the 10 zcro-point-error-added 
AK maps are kept in the final catalog; the number of simulated maps a cluster is found in (out of 
a maximum of ten) is provided in the final cluster catalog, and can be used to select increasingly 
certain subsamples. We use seven detections among the Monte-Carlo maps as our limit because 
the fraction of clusters detected in fewer maps is constant, whereas it drops steeply when more 
detections are required. 

The requirement of > 7 detections from the Monte- Carlo maps removes < 10% of the cluster 
candidates. Because the detection of candidates that appear in fewer maps is very sensitive to 
small photometric errors, these candidates arc likely to be cither false, extremely poor, or at high 
redshift. For the latter two cases, our detection efficiency is very low (as described below), and such 
clusters should not be used in statistical studies. Thus, we expect that this requirement imposes 
no significant bias on the final catalog, and merely serves to reduce the catalog inhomogeneity at 
the limits of our detection procedure. 

5. Contamination, Completeness and Optimizing Detection 

The weakest aspect of optical imaging surveys for galaxy clusters is the lack of distance in- 
formation. The two-dimensional projection of the galaxy distribution can produce many apparent 
overdensities which may be identified as clusters but are not actual physical associations. This 
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issue has been discussed in some detail for the Abeh catalog (see the overview in §2), where rates 
of false clusters due to projection effects as high as 40% and as low as 10% have been claimed. 
Other types of surveys (e.g., X-ray, SZ) arc not as strongly affected by projection, but suffer from 
other selection effects or observational difficulties. Therefore, redshifts are needed to validate the 
existence of a given cluster. 

The number of cluster candidates varies significantly with the detection threshold (Kim et al. 
2002). Obviously, this dependence is also related to the number of false positive detections. As we 
try to achieve high completeness we may suffer from high contamination. Our goal is to minimize 
Nfalse/^candi while maximizing Ncandi where Nf^ise is the number of false positives and N^and is 
the number of candidates for each field. We make use of a simulated background distribution in 
order to evaluate N false, while N^and is derived from the real galaxy catalogs. Because each plate 
suffers from slightly different systematic errors in both photometry and star/galaxy separation, 
our procedure to optimize detection and evaluate completeness and contamination is done for each 
plate individually. This is preferable to selecting a single region of the sky, optimizing the algorithm 
for that region, evaluating the contamination rate, and assuming that is true for the rest of the sky. 



5.1. The Rayleigh-Levy Distribution and Contamination 

Initially, we attempted to estimate contamination rates by producing random catalogs with 
the same mean density as our survey fields. This resulted in obvious underestimation of the 
contamination rate. As pointed out by Postman, Lauer, Oegerle, & Donahue (2002, hereafter 
PLOD), the basic sources for false positives in cluster catalogs based only on photometric data 
are random fluctuations and superpositions of poor galaxy groups. For this survey, in which we 
effectively use a single photometric band for cluster detection, the susceptibility to false detections 
due to projection effects is high. Therefore, we adopt a methodology similar to that of PLOD, 
using the Rayleigh-Levy (RL) distribution to generate simulated galaxy positions. The RL random 
walk creates a linear ordering among the distributed galaxies, implying that galaxies interact only 
with the ones in their immediate vicinity, and each small set of neighbors interacts independently 
of any other small set. The resulting galaxy distribution, as shown by Mandelbrot (1975), follows 
the same two- and three-point correlation functions as described by Peebles & Groth (1975). 

For a given separation 6, the RL distribution gives the probability of finding a pair with a 
larger separation: 

p{>e) = {e/eQ)-'^ if e>eo 

and 

P{>9) = 1 if 9<eo 

We ran a suite of RL simulations with different values for the parameters Oq and d, and chose 
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values for these such that the resulting pair distribution matches the real galaxy distribution for 
the DPOSS data. We use 6q = 500" and d = 0.6. This value of a-lso corresponds to the diameter 
of the smoothing kernel; angular correlations arc thus removed at smaller scales. We find that the 
values of Oq and d do not vary significantly among plates, and therefore adopt the values above for 
all plates. 

For each plate, We create the simulated RL distribution by choosing a random position within 
the plate limits that does not fall into a bad area, and then selecting a separation from the above 
distribution. Prom this new point in the frame we repeat the same procedure until seven galaxies 
are selected. We then choose a new random location and select another new point in the frame, 
repeating the same procedure until seven galaxies are selected. This process goes on until we 
generate a number of positions equal to the total number of galaxies for that particular plate. The 
resulting distributions arc very insensitive to the number of galaxies we choose before restarting. 

We also tested the full random walk distribution, where we start from one location, and 
continue relocating to new points with step sizes following the P(> 6) distribution (as opposed 
to choosing a new random point, distributing a fixed number of points around it, and repeating). 
Again, the final distribution is very similar to the restricted random walk; this is because the 
clustering scale we are interested in is well below the cutoff in the two-point correlation function 
introduced by the restricted random walk. By comparing the angular separation distribution for 
these different sets we find them to be indistinguishable up to scales of 2°. We therefore use the 
restricted RL in our simulations, as it is significantly less computationally intensive than the full 
random walk. 

For each plate in the survey, we produce an RL distribution with the same number of objects as 
the galaxy catalog for that plate, with the same bad areas excised. A density map is generated from 
this distribution, and SExtr actor is run with a range of threshold/minimum area pairs on this map, 
along with the map generated from the real galaxy catalog. We use detection thresholds ranging 
from 900 to 2200 galaxies per square degree, in increments of 20, and minimum areas of 30, 40, 
50 and 100 square pixels (or arcmin). Initial tests using a broader range of parameters indicated 
that the optimal pair would always be found within these ranges. Once detection is completed 
on both the RL and real maps, we measure the contamination rate as C = Nfaise/^cand-, and 
select the parameter set for which C is closest to 10%, while maximizing Ncand'i typically, this 
results in 8% < C < 12%. We performed tests using other values for the contamination rate; 
we found that completeness is not substantially improved by allowing higher C, while lower rates, 
such as C = 5%, produced more widely varying parameter sets due to small fluctuations in C. 
Importantly, enforcing a fixed contamination rate does not produce large plate-to-plate variations 
in the completeness; this can be seen by comparing the completeness functions for a large set of 
plates. Nevertheless, it would be inappropriate to use a mean completeness function for the entire 
survey. 

We then examined the optimal parameters derived for each plate. We found that minimum 
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areas of 30 or 40 square arcminutes produced very similar results, in terms of the distribution 
of both the detection thresholds and final number of candidates per plate. We also examined 
the completeness functions (described below), and found that they were extremely similar for both 
values of the minimum area. Using a minimum area of 40 square arcminutes produced slightly lower 
variance in the number of candidates per plate; we have therefore used this minimum detection 
area for all plates. Physically, this is sensible: we know that the depth and classification accuracy 
vary from plate to plate, which will affect the threshold in galaxies per square degree, but there is 
no reason to expect the sizes of clusters to vary from plate to plate. The larger minimum areas (50 
and 100 square arcminutes) produce radically fewer cluster candidates. 

We stress that the entire optimization procedure is performed separately for each field. The 
optimal parameters for a given field will not be appropriate for another. This is not due only to 
zero point differences or plate variations, but also to the large scale structure (LSS) in the universe 
(Bramel, Nichol &; Pope 2000). Our final cluster catalog is designed to maintain a constant contam- 
ination level across the whole sky. One must also note that the RL distribution does not take into 
account inhomogeneities at large scales; this test is specifically designed to estimate the contamina- 
tion rate due to chance projections of small groups or physically unassociated galaxies. Because it is 
impossible to know a priori which overdensities are true physical associations contributing to large 
scale structure, there is no way to use the actual catalogs to measure contamination rates. Only 
complete spectroscopic follow-up can determine the accuracy of our estimates; as described below 
in §5, our own spectroscopy suggests low contamination rates, consistent with the 10% threshold 
imposed. Additionally, the results of PLOD demonstrate that the RL estimates of false positive 
rates is verified by such spectroscopy. 

An alternative method to estimate contamination rates is to shuffle the galaxy positions in a 

plate catalog, while maintaining the object magnitudes. This method preserves some very large- 
scale correlations, while removing those at smaller scales. As shown in Lopes et al. (2003), this 
technique produces contamination rates (and the corresponding optimized detection parameters) 
which are very similar to those derived using the RL distribution. 

5.2. Completeness From Simulations 

Perhaps the most important feature of this cluster catalog is the detailed assessment of the 
selection function. In order to properly assess this, it is fundamental that each plate be treated 
individually. Because the detection parameters for each DPOSS field are determined separately, 
we assess the completeness as a function of redshift and richness individually for each plate, and 
provide this information for each plate. When generating mock catalogs from simulations (or from 
other observations), one can then apply these selection functions to mimic potential systematic 
effects in our catalog. 

To assess the selection function for each plate, we use simulated galaxy clusters which are 
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placed in the background field derived from the respective plate. We then use these galaxy catalogs 

consisting of a background plus artificial clusters to generate density maps on which we run SEx- 
tractor. When running SExtractor, we use the same detection parameters obtained in the previous 
section, which arc also used for the real cluster detections. 

Artificial clusters are typically defined using a Schccter luminosity function with its choice of a 
and M*; a surface density profile; a given cluster composition; a maximum radius and a core radius. 
Lobe et al. (2000) tested the variation of the completeness rate as a function of the spatial profile 
slope, while PLOD also tested different cluster compositions. As expected, steeper profiles provide 
the highest completeness rates. We discuss the effects of varying different cluster parameters at the 
end of this section. 

An important issue is the background into which we insert the artificial clusters. While the RL 
distribution is well suited to estimate the contamination rate and optimize the detection parameters, 
it only represents an ideal case, being unable to reproduce all the inhoniogcncitics due to large scale 
structure. As we discuss below, the RL distribution overestimates the completeness rate. Thus, 
for each plate, we use the actual galaxy catalog from that plate to provide the background field in 
which we insert the artificial clusters 

The simulated clusters follow a Schecter luminosity function with parameters given by Paolillo 
et al. (2001) (a = -1.1 and M* = -21.53). The constituent galaxies have -23.4 < < -16.4, 
such that our observed magnitude range is fully covered at all redshifts of interest. The surface 
profile adopted is a power law in radius of the form r^, where (5 = —1.3. This is in the middle 
of the observed range, —1.6 < /3 < —1.0, as discussed in Tyson & Fischer (1995) and Squires et 
al. (1996). The galaxies are placed within a maximum radius of Vmax = l-5/i~^ Mpc with a core 
radius of Vcore = 0.15h~^ Mpc. The clusters are composed of 60% elliptical galaxies and 40% Sbc 
galaxies, with SEDs taken from Coleman, Wu, &; Weedman (1980), convolved with the DPOSS r 
filter. Each galaxy has the appropriate ^-correction applied, as well as a random photometric error 
from Gal et al. (2003) added. 

We proceed as follows: 

1. We generate clusters with six different richnesses, which are shifted to seven redshifts, with the 
catalogs trimmed at 15.0 < m,. < 19.5. The richness classes are Ngais = [15, 25, 35, 55, 80, 120] 
galaxies, while the redshifts adopted are z = [0.08,0.12,0.16,0.20,0.24,0.28,0.32]. 

2. For each of the 42 richness/redshift combinations we generate 5 simulated clusters, which are 
placed at random positions in the real background galaxy distribution, avoiding both excised 
areas and cluster candidates detected in the real fields. This is repeated ten times, resulting in 
ten catalogs containing a total of 50 simulated clusters for each richness/redshift combination. 
Thus, there are 420 total simulated catalogs generated for each plate (6 richnesses x 7 redshifts 
X 10 realizations). An example is shown in Figure 2, where we plot Plate 389, with the bad 
areas marked, and five simulated clusters at z = 0.16 and Ngau = 80 inserted. 
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3. We then use the AK to produce density maps for each of these 420 background plus cluster 
galaxy catalogs. This is performed individually for each plate, and SExtractor run with the 
optimal parameters determined from the contamination tests. 

4. Next, we compare the candidate positions in these maps with the initial input positions. We 
use a matching radius of A00h~^ kpc, which is small compared to the typical size of a cluster. 

5.2.1. RL vs. Real Backgrounds 

For one DPOSS field (389) we repeated the entire procedure twice, first using the real back- 
ground and then the RL background distribution. The purpose of this test is to check how well the 
RL distribution reflects the real background. We show the results of this test in Figure 3, where 
the completeness rates obtained using the RL distribution are the dashed lines, while the solid 
lines use the real background. Each panel shows a different richness class, with richness increasing 
from bottom to top. We see that the dashed lines (RL background) are always higher than the 
solid lines, for all richness classes and at all redshifts. At z = 0.16 (the median redshift of the 
DPOSS sample) and Ngais = 35 (close to the median richness value), the recovery rate is 94% using 
the RL background, but only 72% using the real background. This suggests that using the RL 
distribution likely overestimates the completeness. An interesting comparison can be made with 
Figure 11. For each richness class shown there, the median redshift is similar to the redshift where 
the completeness estimate for that subsample starts to decrease significantly. 

We note that Kim et al. (2002) modified their detection limits after comparing the recovery 
rate of artificial clusters using uniform and real backgrounds. Their goal was to improve the 
completeness rate; however, this necessarily has an adverse impact on the contamination. We do 
not apply such corrections, as our goal is to maintain a constant contamination rate. 

5.2.2. Dependence on Cluster Model 

We also wish to test the variation of completeness with the cluster composition, surface density 
profile slope, luminosity function slope, and maximum and core radii. We use a single plate, varying 
the properties of the simulated clusters and repeating our completeness tests. Figure 5 shows the 
dependence of the completeness on the cluster composition and spatial profile slope (/?) . Each row 
represents a different value for the slope (/? = -1.0, -1.3, -1.6, from bottom to top), while each 
column shows the results for different compositions (from right to left: 100% E, 60% E + 40% Sbc, 
and 20% E + 80% Sbc). The improvement when we adopt steeper profiles is evident; as expected, 
tighter cores are easier to detect. There is little if any correlation with the cluster composition. 
This is expected as the difference in the ^-corrections used (for E and Sbc galaxies) is only 0.2"* 
at the highest redshift of the simulations {z = 0.32). For instance, for Ng^is = 35 and z = 0.16, 
using 60% E + 40% Sbc, the recovery rates for /? = (-1.0,-1.3,-1.6) are (62%,72%,80%), respectively. 
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However, when we fix /3 = —1.3, Ngau = 35 and z = 0.16, the recovery rates are (74%, 72%, 70%) 
for cluster compositions of (100% E, 60% E + 40% Sbc, 20% E + 80% Sbc). 

The effects of changing the luminosity function slope, core radius and maximum radius are 
shown in Figure 6. Each row shows the results for one of these 3 parameters. For the richest clusters, 
steepening the LF slope increases completeness, while poor clusters show an opposite trend, albeit 
with large scatter. The variation with core radius shows a clear trend to higher completeness for 
smaller values of this parameter; this is commensurate with the results of the radial profile test. 
Finally, tests with the cutoff radius show the worst recovery rate for r^ax = l-Oh"^ Mpc, while the 
results are about the same for the other two values tested. 

We conclude that our canonical cluster model (cc = — 1.1,M* = —21.53,/? = — 1.3, r^aa; = 
1.5^~^ Mpc, rcore = 0.15^""^ Mpc, 60%E + 40%S) provides a realistic estimate of the catalog 
completeness. For each candidate, we provide a table of the completeness as a function of richness 
and redshift. An example is shown in Table 2 and Figure 4; a similar table is provided online (at 
http://dposs.caltech.edu/dataproducts/) for each plate used in this survey. Additionally, a 
list of the areas excised due to bright stars, airplanes, etc. is provided for each plate at the same 
location. 



5.3. Spectroscopic Confirmation 

To assess the validity of the above results, as well as to calibrate our photometric redshift 
estimator (see the next section and Paper I) and examine the redshift distribution of our clusters, we 
undertook a complete spectroscopic survey of candidates in two DPOSS fields: 447 (14^*30™, +30°) 
and 475 (01^, +25°). These fields were chosen because they are at relatively high galactic latitude 
(+67° and —40°) where the effects of dust are expected to be small, and because scans in all 
three bands were available when this project was started. However, these maps were generated, 
and cluster detection performed, using the methodology discussed in Paper I, which (as discussed 
earlier) differs significantly from the final technique used here. The original list of cluster candidates 
in these two fields was drawn from the individually calibrated plates, using galaxies to rrir = 20. 
Therefore, the sample with follow-up spectroscopy is only a subset of the final cluster catalog in 
these areas, and also contains targets which are no longer candidates. We remind the reader that 
the catalog for Field 475, which is in the SGP, is not a part of the NGP, and thus not in the area 
covered by the catalog presented in this paper. 



5.3.1. Observations 

Spectra of the cluster candidates were obtained between April 1996 and April 2000 with the 
COSMIC instrument (Kells et al. 1998) in re-imaged mode, mounted at the prime focus of the Hale 
5m Telescope at Palomar Observatory. Targets for each mask were chosen from the DPOSS catalog 
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for the appropriate field, selecting only galaxies with rrir < 20 to maintain reasonable exposure times. 
No color selection was applied; this was done to maximize the number of possible objects which 
could have slits placed on them. Multi-slit masks, made of photographic film negatives, with 1.5" 
slit widths, were mounted at the focal plane. The re-imaged pixel scale is 0.399"/pixel; the slits 
therefore correspond to ~3.75 pixels. The available field for spectroscopy is theoretically 8'xl2', 
with the long dimension corresponding to the spatial axis. The field is further limited by the lack 
of an accurate distortion map in the slitmask production software; this requires that slits be placed 
within ~3' of the central axis in the direction perpendicular to the slits. Nevertheless, the 8' x 3' 
field of the masks is larger than the core radius of clusters at all but the lowest redshifts in our 
sample. 

Slitmasks were designed using the cosmicslitmask program, written by M. Pahre. This soft- 
ware simply takes an input list of coordinates, and allows the user to interactively select objects 
for slit assignment. There is no automated slit optimization algorithm, so objects were selected 
visually for each mask, to maximize the number of targets observed per mask. This often resulted 
in densely packed slitmasks, with as many as 45 slits on a single mask, some as short as 10". 

Spectra were taken using a grism with 300 1/mm, blazed at 5500 A, which produced a dispersion 
of 3.O3A pixel"^. The spectra typically cover a wavelength range ~ 3000A to ~ 9000 A, with a 
central wavelength of ~ 6OOOA. Since the dispersing element is a grism and cannot be tilted, the 
wavelength range is fixed and varies only as a result of the vertical displacement of the slits from 
the mask center. Therefore, nearly all spectra are observed from blueward of the 4000A break 
(for ^; = 0), through the redshifted Na lines (5890A rest frame, 7950A at our highest expected 
z = 0.35). The COSMIC CCD has a gain of SAe'/DN and a relatively high read noise of 13e-. 
Exposure times varied greatly from night to night, depending largely on the seeing (which ranged 
from 0.8" to as high as 2.0") and spectrograph focus changes (due to temperature fluctuations). 
Exposures ranged from a single 1800s exposure up to two 3600s exposures per slitmask, with most 
consisting of two 1800s exposures. 



5.3.2. Data Reduction 

All data reduction was performed using the IRAF package. The object spectra were overscan 
subtracted, and flattened using spectra of a halogen lamp reflected off the dome interior. Cosmic 
rays were removed using the szap task (written by M. Dickinson), and night sky emission removed 
by subtraction in the spatial direction on the two-dimensional images, using the background task 
to perform median filtering and sigma clipping. Because the slits were typically very narrow, we 
used only a second order fit along the spatial direction. After sky subtraction, each spectrum was 
traced and optimally extracted in the individual exposures. The extracted spectra were wavelength 
calibrated using spectra of an arc lamp, taken immediately before or after the object spectra, to 
minimize the effects of instrument flexure. The individual extracted, wavelength-calibrated spectra 
were then combined to produce the final spectrum for each object, used to measure the redshifts. All 
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redshifts were measured by RRG using the IRAF task redshift, written by T. Small. This program 
allows the user to visually mark specific spectroscopic features, and input their rest wavelengths. 
The object's redshift is then calculated, and the positions of other common absorption and emission 
lines overplotted. In this way, the user can check the redshift assignments. Such visual measurement 
was necessitated by the poor S/N of the majority of our spectra; automated techniques (such 
as cross-correlations) were strongly affected by residuals from poor sky subtraction. To enable 
completion of the survey in a reasonable time, a trade-off in exposure times vs. number of masks 
observed was made, since we were only interested in measuring redshifts and not any specific 
galaxy properties (such as velocity dispersions or line strengths). Comparison of redshifts obtained 
by different reducers and at different times suggests the redshifts are accurate to ~A2; = 0.002, or 
approximately 600 km s~^. 

In addition, over the course of the survey, our list of candidate galaxy clusters changed as the 
photometric calibration and sample selection improved. This resulted in a significant number of 
targets being observed which are not included in the final cluster sample. This extraneous data 
is actually useful in assessing the rate at which we miss real clusters and is discussed later in this 
section. A total of 3249 individual spectra were visually inspected; of these, 1655 (51%) were 
identified as galaxies, and 326 (10%) as stars. Thus, a total of 1981 spectra were identifiable, a 
61% success rate. Of the galaxies, 1245 (75%) had securely measured redshifts; the remainder had 
insufficient S/N. We also note that the 10% stellar contamination is consistent with the results 
of Odewahn et al. (2003), and should be taken as an upper limit on misclassification, as extended 
sources {i.e., galaxies) tend to have lower S/N spectra, making identification more difficult. This 
large spectroscopic sample constitutes a significant survey in its own right, with other applications 
in addition to those discussed here. 

In Field 447, there arc 64 cluster candidates. Of these, only 24 have usable spectroscopy, with 
all 24 showing evidence for clusters. There are an additional 4 masks targeted at areas with no 
current candidate; these nevertheless show spectroscopic evidence for physical galaxy associations. 
However, they all have z > 0.22 and Ngau < 30, a regime where our completeness is typically only 
~ 30%. There are unfortunately 40 current candidates in this field without spectroscopy. 

The situation for Field 475 is much better, with 55 candidates in our current sample, of which 
37 have corresponding slitmasks. These all show evidence for physical clustering (as discussed 
below). There are 16 masks not associated with any new candidates. Of these sixteen, four are 
near the edge of the plate or the area used for the densitometry spots. There remain twelve masks 
which show evidence for real galaxy associations. Nevertheless, these twelve areas appear to be 
only moderate overdensities; even lowering our detection threshold to allow 30% contamination 
recovers only three of these. Running our richness estimator on these areas, using the spectroscopic 
redshifts, shows that indeed these possible clusters all have Ngau < 30. Additionally, four of these 
12 masks have clusters at z > 0.22. Finally, eighteen current candidates have no spectroscopy; 
these are all of similar overdensity to those candidates with spectra, and therefore likely to be real. 
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The results of the spectroscopic survey are summarized in Table 3. We give the total number of 
candidates, the number (and fraction) with and without spectroscopy, and the number of extraneous 
masks. Details for each candidate with spectroscopic observations arc provided in Tables 4 (for 
Field 447) and 5 (for Field 475). Column (1) provides the candidate name, columns (2) and (3) 
provide the spectroscopic redshift and the number of galaxies at that redshift (if spectra were 
obtained), and column (4) provides comments, if any. 

The above results are shown in Figure 7, for Field 447 (left panel) and 475 (right panel). We 
show the AK galaxy density maps, with blue circles marking the locations of current candidate 
clusters. These clearly correspond to areas with the highest galaxy density. Green circles show the 
locations of slitmasks associated with current candidates, while red circles mark locations where 
spectra were taken but there is no current candidate. 

We show example spectroscopic redshift histograms for some cluster candidates in Figure 8. 
From these redshift distributions, we can clearly sec strong clustering in redshift space, providing 
evidence for true galaxy clusters in most of the slitmasks. The median redshift of our spectro- 
scopically confirmed cluster candidates is Zmed^^-'^- This value is in keeping with the photometric 
redshifts presented in the following section, and with our earlier estimates of the depth of our galaxy 
catalogs. 

5.3.3. Are Redshift Peaks Real Clusters? 

Although the peaks in the redshift distributions provide evidence for clustering, we would 
like to quantify the likelihood that a given peak is a real structure. Certainly, large peaks (as in 
candidate NSC005559+262442, with 24 galaxies) are undoubtedly clusters. However, many of the 
masks show peaks with only 3 or 4 members. Individual galaxies are assigned as members of a 
peak if they fall within Az = 0.005 (±1500 km s"-*^) of the peak center. Increasing this range does 
not change the results significantly, since most redshift peaks are well isolated. 

A rough estimate of the likelihood for such small peaks being real can be made following the 
argument of Zaritsky et al. (1997). If the underlying redshift distribution of galaxies were smooth, 
the formal probability of a pair of galaxies within Az = 0.005 would be negligible, much less for 
three or more galaxies. However, because galaxies are correlated, the likelihood of pairs, triplets, 
etc., is increased. If we assume that all slitmasks where the richest clump has only three or four 
members correspond to spurious detections (there are 18 such masks out of 95 total masks), the 
probability of finding a single false peak is 18/95, or ~20%. We might then expect a comparable 
fraction of our histograms to show a second redshift peak with similarly few members, but there 
are only eight (~8%). In fact, some of the distributions with multiple peaks may be projections of 
multiple poorer clusters which we detect as a single candidate. Additionally, there is only one mask 
with a secondary peak containing five members. An additional argument in favor of the reality of 
three-or-more galaxy clumps is that if peaks with three members were often spurious, we expect 
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to find many more such peaks than those containing four galaxies. However, we find an almost 
identical number of peaks with three and four galaxies, supporting the idea that even thrcc-galaxy 
peaks are often real. Therefore, we conservatively identify all peaks with four or more members 
as real clusters, and peaks with only three members as tentative. The number of clumps with 
two galaxies is much higher (nearly every mask contains one); we therefore ignore those completely. 
Furthermore, while Zaritsky et al. follow similar arguments for only two galaxies, their spectroscopy 
covers a significantly larger magnitude range (m^ < 22 compared to our rrir < 20) , thereby greatly 
increasing their likely contamination rate. Conservatively, we estimate that at least 80% of our 
robust candidates arc true clusters; this is roughly consistent with our attempt to generate catalogs 
with ~ 10% contamination. 



6. Photometric Redshifts and Richnesses 

Two of the most fundamental properties of galaxy clusters are their distance (or redshift), and 
mass. Because the latter is difficult to measure accurately, and impossible from our photometric 
data alone, we instead measure the richness. The details of our redshift and richness estimation 
techniques are given below. It is important to note that both are measured in a completely model 
independent way; this is possible because, these properties are measured after detection is com- 
pleted, and arc not an inherent part of the procedure. Wc note that this is inherently different 
from many other techniques; for instance, the Matched Filter simultaneously performs detection 
and redshift/richness estimation, while others, such as the Cut Sz Enhance method (Goto et al. 
2002) or the RCS (Gladders &; Yee 2000) output the redshift, based on galaxy colors used as part 
of the detection procedure. 



6.1. Redshift Estimates 

Paper I presented an extremely simple yet effective photometric redshift estimator for DPOSS 
cluster candidates, based on the assumptions that each candidate is a single cluster, at one red- 
shift, and the cluster galaxy population is dominated by early-type galaxies. The estimator is an 
empirical relation between spectroscopic redshift and the median g—r color and mean r magnitude 
of the background-corrected galaxy population for each candidate. We use both colors and magni- 
tudes as we find both to be equally well correlated with the spectroscopic redshift. We count the 
number of galaxies as a function of color, Ng^r, and the number as a function of r magnitude, A^^; 
inside a radius of 1 Mpc (0.67x i2^f,eii)- The background galaxy color and magnitude distributions 
{^bg,g-r ^nd Njjg^^) are determined independently for each plate, scaled to the appropriate area, 
and subtracted from the color and magnitude distributions of each candidate cluster. The median 
g—r color and mean r magnitude of the remaining galaxies is then calculated. This differs slightly 
from the methodology in Paper I, where universal background distributions, taken from a large 
contiguous area, were utilized. In practice, this must be an iterative procedure, because we do not 
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initially know the redshift, and therefore the Abell radius, for our cluster candidates. We therefore 

start with an initial radius corresponding to 1 Mpc at z = 0.05, within which we measure the 
above quantities and derive the initial redshift estimate. The Abell radius is adjusted using this 
new redshift estimate, and the procedure repeated until the redshift converges. 

Spectroscopic redshifts for the clusters used to derive our empirical photometric redshift esti- 
mator were taken from Struble & Rood (1999, hereafter SR99). We do not use our own spectroscopic 
data, as they cover only two plates, and we do not want to bias our estimator. The SR99 data 
are culled from diverse sources in the literature, and are of highly variable quality. A majority 
of the cluster redshifts, especially those at 2; > 0.1, are based on very few (one or two) individual 
galaxy redshifts, and may therefore be incorrect. Nevertheless, this catalog is the largest, somewhat 
homogenized resource currently available. The larger area utilized in this paper results in a final 
training sample comprising 369 clusters (compared to only 46 in Paper I), with a median redshift of 
Zmed = 0.138. Because we now have many more clusters, and their redshift distribution is heavily 
weighted with low-redshift clusters, we bin the colors and magnitudes into spectroscopic redshift 
bins. We divided the range 0.02 < 

Zspec ^ 0.35 into ten bins, and calculate the median Zspecj 
{9 ~ f)medj a-nd rmean for each bin. These binned values are used to derive an empirical relation 
between redshift, median g—r color, and mean r magnitude using a bivariate least-squares fit: 

zphot = 0.5694 X (5 - r)med - 0.00215 x r„ - 0.057243 (1) 

The upper panel of Figure 9 shows the photometric redshift against the spectroscopically mea- 
sured redshift for these 369 Abell clusters, with the lower panel showing the residual. Performing 
the complete iterative procedure, where we begin without assuming the spectroscopic redshift to 
define a starting radius, yields a Qsigma of (zspect - Zphot)/i^ + Zspec) = Az = 0.033. 

In addition, we tested the effect of varying the starting cluster radius from 5' (corresponding 
to ~ l.OMpc at 2; = 0.2) to 15' (~ l.OMpc at z = 0.05). In principle this could lead to different final 

photometric redshifts if the convergence procedure is unstable. The results are shown in Figure 10; 
the Qsigma for Azphot = 0.004. There are a very few outliers, which demonstrates the robustness of 
our simple technique. This also suggests that estimating Zp^ot with two different initial radii and 
selecting those with varying results can be used to test for clusters where the estimated redshift is 
questionable. The starting redshift also does not affect the accuracy of the estimator. 

The photometric redshift estimates for the candidate clusters are provided in the fourth column 
of Table 6. Those clusters where the photometric redshifts using the two starting redshifts disagree 
by more than s/O, times the redshift error are marked with a colon; in cases where the estimator 
failed entirely, we set z = 0.0. The mean redshift for our sample is Zphot,med = 0.1579, which is 
comparable to our original estimate of <z>~0.15, based on the magnitude range covered. Our 
highest redshift clusters are at z = 0.3, with many at z> 0.2. The redshift distribution of our robust 
sample is shown in Figure 11. The entire sample is shown in the top panel, with subsamples of 
varying richnesses in the lower bins. The median redshift for each sample is marked. As expected, 
poor clusters are found only at low redshift, while richer clusters can be seen to larger distances. 
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6.2. 



Richness 



A fundamental physical property of galaxy clusters is their mass. Mass can be measured using 
X-ray data, velocity dispersions, lensing (strong and weak); all of these methods require data which 
is difficult to obtain for extremely large cluster samples, and often impossible for low-mass systems. 
One expects that the number of galaxies in a chistcr (the cluster richness) should be correlated 
with the total cluster mass. Unfortunately, cluster richnesses are notoriously difficult to measure, 
and there are a variety of different estimators. Abell's richness, for instance, is extremely poorly 
correlated with mass; other estimators fare somewhat better (Yee &; Lopez-Cruz 1999). Because 
our catalog covers a large fraction of the sky, we can measure optical richnesses for many clusters 
which have spectroscopic data, as well as prepare specific subsamples for future spectroscopic study. 

Given a redshift estimate for each cluster, we can measure a richness that samples the same 
absolute magnitude range. We compute the richness in a fixed absolute magnitude interval, M*—l< 
M <M* + 2, using the r-band data, and assuming M* = —21.53, taken from the cluster luminosity 
functions derived by (Paolillo et al. 2001). Galaxies are counted within a radius of lh~^ Mpc. For 
clusters at low redshift (z < 0.17), we measure the richness by directly summing the background- 
corrected number of galaxies in the appropriate magnitude interval. Only at these low redshifts 
is this entire magnitude range contained within our data. For more distant clusters {z > 0.17), 
we directly sum the galaxies with M* — 1 < M < M20, where M20 is the absolute magnitude limit 
corresponding to our catalog limit = 20.0. We then calculate a correction factor 7 for the 
richness, defined as 



We set the exponent in the luminosity function to a = —1.1, as we do for our simulated clusters; 
the correction factor is slightly larger for lower a (by ~5 — 10% for a = —0.87). Typical values of 
7 are 1.125 for z = 0.2 and 1.50 at z = O.S. 

We use a background determined individually for each plate. We simply take the distribution of 
galaxies as a function of magnitude for the entire plate, and subtract this (scaled to the appropriate 
area) from the same distribution for each cluster. This ensures that any systematic errors in 
an individual plate (calibration, classification) are maintained when performing the background 
subtraction. Our methodology lies between completely local background estimation (such as taking 
an annulus near the cluster), and completely global estimates (which would use the entire survey 
area) . Tests using median filtering to remove overdense regions from the background measurement 
show that this method does not change the final richnesses appreciably; the mean richness is 
increased by ANgais = 3) if this alternative method is used. 

The final richness {Ngais) for each cluster is listed in the fifth column of Table 6. We show 
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the richness distribution of our robust cluster candidates in Figure 12, for the whole sample and 
in different redshift bins. At high z, where our completeness drops, we find only richer clusters, 
while at low z, where the volume is smaller, we find mostly poor clusters. The median richness for 
our clusters is Ngais,med = 31, which corresponds to the poorest end of Abell's richness class 0. We 
therefore expect that our catalog contains many very poor clusters (or even groups). 

6.2.1. Richness Errors 

Errors in our photometric redshift estimator will introduce an error in the measured richness 
for each cluster. We investigated the magnitude of this effect by calculating richnesses for the 
Abell clusters with spectroscopic rcdshifts. We calculate the richness using both the spectroscopic 
redshift and our photometric redshift. The results are shown in the upper left panel of Figure 13. 
The two estimates are very well correlated, with no systematic offset, and larger scatter at lower 
richness, where the redshift estimates are likely to be worse. 

We also tested various cluster radii for measuring richness. The upper right panel of Figure 13 
shows the ratio of richnesses using 1.5Mpc and IMpc radii, versus the l.SMpc radius richness. We 
see that the ratio is approximately constant regardless of richness, and that the ratio is ~ 1.3, less 
than the ratio of the areas. This is clearly because the outer parts of the clusters have low density. 
We therefore use the IMpc radius richnesses, as the background subtraction over a larger area will 
introduce more noise. 

There are many other potential sources of error, both random and systematic. A future paper 
(Gal & de Carvalho 2003) will address many of these issues in detail, including the choice of 
background area, cosmology, fc-corrections, the radius used, etc. We will measure cluster richnesses 
using a variety of techniques; this will allow detailed tests to find which methodology provides the 
best surrogate for a mass measurement. Nevertheless, we are confident that our ad hoc estimator 
provides a useful measurement. We have compared our richness estimator to the A^; measurement 
of Kim (2001) for clusters detected by both surveys. The results are shown in Figure 14, where solid 
circles are clusters where the redshift estimates agree (within errors), while open circles are clusters 
where they do not. The two richness measures are extremely well correlated, suggesting that they 
are both measuring similar properties of the clusters. The agreement is quite remarkable, as they 
are based on different surveys, and measure different properties. One method (Ngais) simply counts 
galaxies, while the other considers the luminosity of the galaxies (Ad). 

6.2.2. Comparison to Abell 

Finally, we attempted to compare our richness estimate to those of Abell. The results are 
shown in the lower left panel of Figure 13, with each comparison cluster as a small dot, and the 
boxes representing binned data. Although a correlation is evident, the scatter is large; significantly 
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larger than the errors due to our measure alone. Abell's richness estimate, counting from the 
third brightest galaxy to two magnitudes fainter, suffers from many drawbacks. First, his redshift 
estimates are often incorrect, causing errors in the radius used to count galaxies. The use of 
an apparent magnitude, ma, to set the magnitude limits means that different clusters have their 
richnesses measured in different absolute magnitude ranges, and introduces random errors due to 
bright foreground galaxies. This effect is shown in the bottom right panel of Figure 13, where 
wc plot Abell's richness against our estimate of the cluster richness measure using a version of 
his technique (detailed in Kim 2001). A similarly poor, but extant, correlation can be seen. The 
difficulties with Abell's richness measures are well documented, as discussed in §1. 

7. The Cluster Catalog 

The first installment of the Northern Sky Optical Cluster Survey catalog presented here covers 
5834n°, and contains 8,155 cluster candidates, yielding 1.4 clusters per square degree. In com- 
parison, Abell's 1958 survey covered a much larger area, yet contains only 2,712 clusters. Our 
catalogue thus represents an increase of roughly one order of magnitude in cluster counts. More 
importantly, these clusters have been selected using an automated, objective algorithm, with ex- 
tremely well characterized contamination rates and selection functions. Each cluster also has a 
consistently measured richness and photometric redshift. Figure 15 shows the sky distribution of 
all our cluster candidates. The median redshift of our sample is z^ed = 0.1579, with a median 
richness of Ngais,med = 31. Thus, our sample is both somewhat deeper than Abell, and extends to 
significantly poorer systems. As mentioned earlier, the redshift distribution is shown in Figure 11, 
and the richness distribution in Figure 12. DPOSS i^-plate images (500" diameter) of four new, rich 
clusters are shown in Figure 16. Two were not previously known, and two (NSC172013+264028 = 
RXC J1720.1+2637 and NSC122906-f-473720 = RXC J1229.0+4737) were detected only as X-ray 
clusters. 

We provide the complete cluster catalog for our survey area in Table 6, as well as at 
http://dposs.caltech.edu/dataproducts/, sorted by increasing RA. 

The table provides: 

1. The cluster name. The naming convention is NSChhmmss-|-ddmmss. Coordinates are J2000. 

2. The cluster RA and 

3. The cluster Dec, in decimal coordinates, also J2000. 

4. The photometric redshift. If this is zero, the estimator failed to converge. 

5. The measured richness. If Zphot > 0.17, the correction factor 7, as described in §6.2, is applied. 
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6. The number of simulated maps in which the cluster was detected. Only clusters detected in 
seven or more maps are included here. 

7. The plate number from which the cluster is drawn. 

7.1. Consistency with Previous Surveys 

Although our data and methodology is different from those of prior surveys, we expect that 
optically-based cluster searches will find similar objects. Thus, we may ask if our results are 
consistent with other recent surveys. 

The obvious survey to compare with is that of Abell. In the area covered in this paper, there are 
1090 Abell clusters. Of these, ~ 93% are recovered by us; an exact number is difficult to determine 
due to centroiding errors, varying matching radii, unknown redshifts, proximity to excised areas, 
and other effects; some of these will be treated in a full comparison when the second (and final) 
portion of our catalog is completed. Nevertheless, the high recovery rate suggests that the Abell 
catalog consists mostly of real clusters (as argued by Miller et al. (1999)), but is incomplete, as our 
sample contains eight times more clusters. Unfortunately, it is difficult to assess the Abell catalog 
incompleteness in the regime where it is used for statistical studies (A^yi6eH > 50, z < 0.2) because 
of the poor correlation between Abell's richness and our Ngais- If we assume Ngais and N Abell are 
equivalent, then the two catalogs contain very similar numbers of clusters in the aforementioned 
range. However, a small shift in the limits used for comparison (for instance, setting Ngau > 45) 
would lead one to conclude that the statistical Abell catalog is ~ 35% incomplete. These discrepant 
results are consistent with the difficulties encountered by other authors in assessing the completeness 
of the Abell catalog. 

An excellent alternative comparison sample is that of PLOD, who discuss the surface density 
of clusters as a function of redshift. They find 1.19 clusters per square degree for z < 0.2 and 
Ac/ > 40. In our 5800 square degree sample, we therefore expect 6902 clusters for the same range 
of parameters. Unfortunately, we do not measure A^;, but PLOD provide a conversion to Abell 
richness, A^,; = 1.24 x RAbell'i their A^,; cut thus corresponds to RAbeli ~ 30. At these low richnesses, 
the Abell richness corresponds roughly to our Ngais (see Figure 13). We therefore apply these 
cuts {zphot < 0.2, Ngais > 30) to our catalog, resulting in a sample of 2888 clusters. We then 
apply our selection function (shown in Figure 3); at the higher end of the redshift range (z ~ 0.15) 
our detection efficiency is quite low for the most common, poorer clusters. The estimated actual 
number of clusters expected in our survey area, after applying this correction, is 6461 clusters, 
compared to 6902 predicted by the constant comoving space density given by PLOD. This is only 
a 6% difference, which is remarkable given that both surveys use different techniques, and have to 
apply estimated selection functions to arrive at actual space densities. 

This excellent agreement suggests that both our survey and the smaller but deeper survey of 
PLOD detect similar objects in the overlapping redshift range. In addition, the selection functions 
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presented by the two surveys provide realistic estimates of their completeness. 

We note also that PLOD find a space density of Abell-like clusters that is a factor of 1.5 higher 
than in the Abell catalog. This result relies on a conversion of to NAbelu which, like our own 

conversion, is somewhat uncertain. Nevertheless, the agreement in density between our survey 
and PLOD does suggest significant incompleteness in the Abell catalog. Because the number of 
detected clusters is strongly increasing at the poor end of Abell's R = 1 bin, these results are very 
sensitive to small changes in the richness cuts applied. 

A preliminary comparison to the NORAS X-ray cluster survey (Bohringcr ct al. 2000) shown 
in Gal (2001) demonstrates that nearly 90% of NORAS clusters aX z < 0.2 are detected in DPOSS. 
In contrast, the DPOSS cluster catalog contains over an order of magnitude more clusters than 
NORAS; in most cases these are poor systems. To understand the biases in these catalogs, we 
are undertaking a full comparison in which the same properties {fx, J^gals, ^phot) consistently 
measured for both NORAS and DPOSS clusters. 

8. Discussion and Future Directions 

Wc have presented the community a large, new cluster catalog, which meets the original 
objectives outlined at the beginning of this paper. This catalog includes accurate photometric 
redshifts and richnesses for a sample of 8,155 cluster candidates. Additionally, we provide detailed 
selection functions in both redshift and richness, on an individual plate basis, which can be used in 
the generation of mock catalogs for testing cosmological models, measuring the correlation function, 
and other large scale structure studies. 

Certainly, this will be superseded in the future by deeper surveys, with more accurate pho- 
tometry in more bandpasses, and with superior spectroscopic or photometric redshift information. 
Our catalog is still limited by the use of the projected galaxy distribution to detect clusters, and 
the poor photometric accuracy in comparison to modern digital surveys, such as the SDSS. 

In the future, we will publish an additional catalog, covering the Southern Galactic Cap region, 
and including the less well calibrated plates from the NGP. The final Northern Sky Cluster catalog 
will then cover over 10,000 square degrees. 

The scientific uses for this catalog are numerous. Measurements of statistical properties of the 
cluster population, such as the cluster-cluster correlation function ^cc and the cluster mass function 
are all avenues of further research. Comparisons to the relatively large X-ray selected samples from 
the RASS are forthcoming, with the potential to find many poorer systems by performing joint 
optical/X-ray detections. The existence of independent samples of thousands of clusters will allow 
us to evaluate the biases present in the different detection methods. 

We also plan to examine the multiplicity function (Puddu et al. 2003), from small groups to 
rich clusters. This requires a separate catalog of poor groups, which is being constructed using 
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a different algorithm by (lovino et al. 2003). The catalog presented here and the compact group 
catalog are expected to have significant overlap, which will also provide important cross-checks in 
the difficult domain of small galaxy associations. 

This sample is also fertile ground for the selection of specific, well-defined subsamples for 
follow-up observational studies. Using these cluster locations and the now public DPOSS data, we 
can search for clusters with substructure, excess blue galaxy populations, unusual optical/X-ray 
flux ratios, or any of a variety of interesting properties. Other methods for cluster detection have 
also been applied to our galaxy catalogs (Puddu et al. 2001), which can be used for further testing 
and comparisons. 

We remind the reader that all of the data for this survey (including future installments) can 
be found at http://dposs.caltech.edu/dataproducts/. 
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Fig. 1. — The effect of varying the initial smoothing window on cluster appearance. Each panel con- 
tains a simulated background with four simulated clusters, as described in the text. The smoothing 
kernel ranges in size from 300" for the upper panel, to 800" for the lower panel, in 100" increments. 
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Fig. 2. — Plate 389, with the bad areas marked, and five simulated clusters at z = 0.16 and 
^gals = 80 inserted. 
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Fig. 3. — The completeness rate evaluated with two different backgrounds: a RL distribution 
(dashed line) and a real galaxy catalog from DPOSS (full line). The richness, increasing from 
bottom to top, is indicated on each panel. 
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Fig. 4. — An example eompleteness function, for Plate 389. From top to bottom, the functions 
correspond to clusters of richnesses Ngau = [120, 80, 55, 35, 25, 15]. 
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Fig. 5. — The selection function evaluated with different cluster compositions and spatial profile 
slopes p. From bottom to top /3 assumes the values -1.0, -1.3 and -1.6. The cluster composition, 
from right to left panels is 20% E + 80% Sbc, 60% E + 40% Sbc, and 100% E. Richness classes 
are the same as Figure 4. 



-35- 



> 
o 
o 



90 
60 
30 

90 
60 



maximum radius = 1 .0 Mpc 



luminosity function slope = -0.80 



maximum radius = 1.5 Mpc 




luminosity function slope = -0.95 



luminosity function slope = -1.1 



O 
CD 

1 30 

(U 
CL 





90 
60 
30 


.1 .2 .3 .1 .2 .3 .1 .2 .3 

redshift 

Fig. 6. — Dependence of the selection function with the luminosity function slope (bottom panels) , 
core radius (middle panels) and cutoff radius (top panels). The values tested, from left to right, 
are a = —0.8,-0.95,-1.1 for the LF slope, rcore = 0.1, 0.15, 0.2/i^^ Mpc for the core radius, and 
i^max = 1.0,1.5,2.0 Mpc for the maximum radius. Richness classes are the same as shown in 
Figure 4. 
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Fig. 7. — AK galaxy density maps for Fields 447 (left) and 475 (right), with blue circles marking 
the locations of current candidate clusters, green circles showing slitmasks associated with current 
candidates, and red circles showing locations where spectra were taken but there is no current 
candidate. 
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Fig. 8.— 



Spectroscopic redshift histograms for selected candidates in Field 475. 
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Fig. 9. — The photometrically estimated redshift vs. the spectroscopically measured redshift for 
369 Abell clusters. Residuals as a function of redshift are shown in the bottom panel. 
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Fig. 10. — The photometrically estimated redshift using starting radii of 5'and 15'. Differences 
between the estimates are shown in the bottom panel. 
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Fig. 11. — The redshift distribution of our candidate clusters. The distribution for all clusters is 
shown at top, with lower panels dividing the sample by Ngais- 
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Fig. 12. — Richness distributions for the entire sample (top panel) and in different redshift bins 
(lower panels). 
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Fig. 13. — Richness estimator tests using Abell clusters with spectroscopic rcdshifts. Panel (a) 
shows our richness measure using the spectroscopic redshift from SR99 compared to the richness 
using our photo-z. Panel (b) is the ratio of richnesses using l.SMpc and IMpc radii. Panel (c) shows 
our richness (using Zgpec) against Abell's. Panel (d) compares Abell's richness to our measurement 
using the Kim (2001) technique to count galaxies with rris < m < ms + 2.1n the lower panels, each 
small dot is a single cluster, with the large squares representing binned data. 
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Fig. 14. — Comparison of our richness measure, Ngais, with independent measurements of from 
Kim (2001). 
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Fig. 15. — The sky distribution of our candidate clusters. Black lines show the boundaries of plates 
used in this paper; dotted boundaries are unused plates with poorer calibration. 




Fig. 16. — DPOSS F-plate images of four rich clusters. Two were not previously known, and two 
(NSC172013+264028 = RXC J1720.1+2637 and NSC122906+473720 = RXC J1229.0+4737) were 
detected only as X-ray clusters. 



-46- 



Table 1. Percentage of Clusters Detected vs. Number of Maps 



Percent 


N(maps) 


Percent 


N(maps) 


99.9 


1 


92.8 


6 


99.7 


2 


89.7 


7 


98.9 


3 


68.6 


8 


98.5 


4 


66.3 


9 


96.3 


5 


62.1 


10 



Table 2. Example Completeness Function: Plate 389 



N/z 


0.08 


0.12 


0.16 


0.20 


0.24 


0.28 


0.32 


15 


0.16 


0.14 


0.08 


0.04 


0.04 


0.02 


0.00 


25 


0.64 


0.46 


0.38 


0.24 


0.00 


0.04 


0.04 


35 


0.94 


0.78 


0.76 


0.28 


0.16 


0.06 


0.00 


55 


0.98 


0.94 


0.84 


0.78 


0.34 


0.08 


0.10 


80 


1.00 


1.00 


1.00 


0.86 


0.64 


0.32 


0.00 


120 


1.00 


0.98 


1.00 


0.94 


0.88 


0.46 


0.38 



Table 3. Spectroscopic Survey Results 



Field 


N(cands) 


N(spec) [frac] 


N(extra) 


447 


64 


24 [0.38] 


4 


475 


55 


37 [0.67] 


16 
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Table 4. Field 447 Spectroscopy 



Candidate/Mask 


^spec 


N{z) 


Comment 


NSC142311+320840 


0.200 


10 




NSC142841+323859 


0.127 


5 




NSC142920+270609 


0.268 


6 




NSC142937+301403 


0.103 


4 




NSC143203+293404 


0.221 


14 




NSC143237+313532 


0.131 


13 




NSC143330+292738 


0.219 


12 




NSC143400+301222 


0.222 


3 




NSC143437+284024 


0.205 


11 




NSC143539+281143 


0.203 


3 




NSC143737+300923 


0.338 


9 




NSC143744+302547 


0.160 


13 




NSC143910+290229 


0.253 


11 




NSC144210+294444 


0.216 


4 




NSC144229+292545 


0.224 


8 




NSC144231+323227 


0.243 


11 




NSC144250+314342 


0.243 


9 




NSC144315+305758 


0.227 


4 




NSC144432+311149 


0.233 


14 




NSC144457+300112 


0.177 


3 




NSC144603+301148 


0.109 


5 




NSC144635+281740 


0.229 


4 




NSC144713+302554 


0.170 


4 




NSC144820+272134 


0.233 


20 




144328+313136 


0.240 


11 


High z\ Ngais = 36.3 


144352+302724 


0.320 


4 


High z; Ngais = 22.5 


144404+313214 


0.233 


6 


High z; Ngais = 25.3 


144902+323713 


0.200 


3 


Ngals = 13.2 
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Table 5. Field 475 Spectroscopy 



Candidate/Mask 




N{z) 


Comment 




1 QA 


94 








o/o 








1 fl /fi 








1 1 

1 1 








7 




N^inni n9'ii -i-9'i9n9s 


. i-OiJ 


1 <^ 






.zuu 


u 






1 fifi 

. luu 


1 9 

i.Z 






1 97-- 


Q 
o 




N^Pni n'^4.S-l-9fi9fi9S 


949 


1 9 
±z 




N^Pfii n4n'^-i-9^^Qnfi 

iNDv^uiu^uon^zooyuu 


. Z4:0 


7 




NSPfii n4ns-i-9^nfi^zL 

iM Dv^U±U^Uon^ZtJUUtJ^ 




1 7 




NClPfll (14904-971 898 




e 
o 




NClPfll n4^4-l-9fi'?Qn8 




7 




iNo^uiu^oyn^zo^uio 


1 ^8 


1 






941 


1 7 






1 fi4 


1 11 




MC;Pm n74S-U9zL'?791 


.z,oo 


7 






1 Q9 


o 




N'^Pm n7^S-l-979fi9fi 


117/ 9'^J^ 

.LL 1 / .ZOO 






N^IPni n84.7-l-9'i991 A 


900 


1 9 
±z 




N^pni nsl^'^-^9zlf^'^1 1 

INuv^UlUOOOT^Z^OOll 


1 Q7 








11*^ 


7 




INJCiPni 1 1 i^9j-974fi1 9 




o 
o 




N^Pfll 1 91 n-l-949fi4.fi 




1 7 




N^pni 1 9^11 -i-9'infini 

IN kJ^Ul IZOln^ZOUOUl 


1 

.loo 






NSPfll 1 444-1-9449 tin 


1 7R 

. -L 1 O 


a 
o 




NCiPfll 1 ^^91 -l-949Q9ii 
IN o v^u ± ±oz in^z'izyztj 




R 

u 




N^Pni 1 fini -l-9997'^fi 

iNuV./UllUUln^ZZZ(OU 


94R 

. Z4:0 


7 
t 




N^Pni 1 fim -l-97'^Q'?8 
iNov./Uiiuuin^z(oyoo 


1 90 
. J_ zu 


Q 
o 




NSCOl 1726-1-225257 


.123 


6 




NSC011825+273800 


.177 


10 




NSC011932+243213 


.142 


7 




NSC011954+244954 


.207 


5 




NSC012049+233053 


.117 


21 




NSC012057+245751 


.190 


9 




010030+243949 


.082 


7 


Ngals = 17.1 


011252+222124 


.140 


6 


iVga/s = 16.0 


010432+243635 


.266 


9 


Ngals = 20.4 


005940+262717 


.193 


14 


Ngals = 22.8 


005645+230911 


.180 


6 


iVga/s = 18.9 


010441+222618 


.250 


5 


High z; Ngals = ■ 
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Table 5 — Continued 



Candidate /Mask 


^spec 


Niz) 


Comment 


010752+253143 


.200 


6 


Ngals = 29.4 


011142+230337 


.195 


7 


Ngal. = 9.8 


011818+224753 


.267 


5 


High z; Ngals = 27.63 


010025+242744 


.125 


10 


Ngals = 20.3 


011000+250756 


.200 


6 


Ngals = 26.7 


011050+231344 


.115 


5 


Ngals = 23.4 


010027+251510 


.227 


5 


High z; Ngals = 28.6 


011548+252229 


.185 


9 


Ngals = 22.8 


011747+245810 


.145 


11 


Ngals = 23.2 


005608+251240 


.175 


8 


Ngals = 16.80 



Table 6. The Northern Sky Cluster Catalog: Excerpt 



Name 


RA (J2000.0) Dec (J2000.0) 


Zphot 


Ngals 


maps 


Plate 


NSC080140+623236 


120.41708 


62.54339 


0.0732 


28.1 


10 


124 


NSC080142+575020 


120.42686 


57.83877 


0.1350 


20.6 


7 


124 


NSC080205+634700 


120.52274 


63.78326 


0.0725 


45.9 


10 


089 


NSC080217+633240 


120.57254 


63.54455 


0.0626 


22.9 


10 


089 


NSC080218+604112 


120.57571 


60.68660 


0.1080 


23.6 


10 


124 


NSC080244+642523 


120.68704 


64.42318 


0.1107 


28.6 


7 


089 


NSC080245+612414 


120.69027 


61.40387 


0.1149: 


29.8 


10 


124 


NSC080300+631419 


120.75366 


63.23854 


0.1434 


38.3 


8 


089 


NSC080320+652141 


120.83540 


65.36151 


0.1470 


29.5 


7 


089 


NSC080327+640531 


120.86494 


64.09192 


0.1341 


28.6 


7 


089 


NSC080337+585755 


120.90653 


58.96515 


0.2260 


23.4 


10 


124 


NSC080347+630320 


120.94792 


63.05546 


0.0874 


30.0 


9 


089 


NSC080401+625402 


121.00716 


62.90060 


0.0850 


33.9 


8 


089 


NSC080435+395527 


121.14706 


39.92425 


0.1014 


27.7 


10 


312 


NSC080542+670818 


121.42574 


67.13827 


0.0942 


32.2 


10 


089 


NSC080608+374711 


121.53357 


37.78649 


0.2323: 


30.2 


10 


312 


NSC080615+381832 


121.56293 


38.30891 


0.2619: 


60.7 


10 


312 


NSC080621+614129 


121.58971 


61.69127 


0.1691 


33.1 


10 


124 


NSC080650+620358 


121.71079 


62.06607 


0.2518 


37.6 


10 


124 


NSC080717+ 393521 


121.82310 


39.58903 


0.0936 


18.6 


10 


312 


NSC'()8()715-(i()r)()2<) 


121.91125 


60.811:57 


0.2797 


29.7 


10 


124 



